Quantum probe and design for a chemical compass with magnetic nanostructures 
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Magnetic fields as weak as Earth's may affect the outcome of certain photochemical reactions 
that go through a radical pair intermediate. When the reaction environment is anisotropic, this 
phenomenon can form the basis of a chemical compass and has been proposed as a mechanism 
for animal magnetoreception. Here, we demonstrate how to optimize the design of a chemical 
compass with a much better directional sensitivity simply by a gradient field, e.g. from a magnetic 
nanostructure. We propose an experimental test of these predictions, and suggest design principles 
for a hybrid metallic-organic chemical compass. In addition to the practical interest in designing a 
biomimetic weak magnetic field sensor, our result shows that gradient fields can server as powerful 
tools to probe spin correlations in radical pair reactions. 
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Introduction. — Recently, there has been increasing in- 
terest in quantum biology namely investigating quantum 
effects in chemical and biological systems, e.g., light har- 
vesting systems [l|, avian compass 2-5] and olfactory 
sense [6|. The main motivation is to understand how 
quantum coherence (entanglement) may be exploited for 
the accomplishment of biological functions. As a key step 
towards this goal, it is desirable to find tools that can 
detect quantum effects under ambient conditions. The 
ultimate goal of practical interest in studying quantum 
biology is to learn from nature and design highly effi- 
cient devices that can mimic biological systems in order 
to complete important tasks, e.g. collecting solar energy 
and detecting weak magnetic held. 

As an example of quantum biology, the radical pair 
mechanism is an intriguing hypothesis Q to explain the 
ability of some species to respond to weak magnetic fields 
T3|, e.g. birds [UHl!, fruit flies [H], and plants 
A magnetochemical compass could find applications 
in remote magnetometry, in magnetic mapping of mi- 
croscopic or topographically complex materials, and in 
imaging through scattering media [l7j |. It was demon- 
strated that a synthetic donor-bridge-acceptor compass 
composed of a linked carotenoid (C), porphyrin (P), and 
fullerene (F) [l8| can work at low temperature (193 K). It 
is surprising that such a triad molecule is the only known 
example that has been experimentally demonstrated to 
be sensitive to the geomagnetic field (yet not at room 
temperature). It is currently not known how one might 
construct a biomimetic or synthetic chemical compass 
that functions at ambient temperature. 

In this Letter, we approach to the goals of studying 
quantum biology in the context of chemical compass by 
demonstrating that a suitably designed gradient field can 
significantly improve the performance of a model chemi- 
cal compass (apart from increasing the intersystem cross- 
ing rate [H}), see Fig.Q] It also opens a possible route to 
probe spin correlations of radical pairs and thereby inves- 
tigate the role of quantum effects in spin chemistry. The 
gradient field is strong at the location of one spin, and ap- 
proximately zero at the other. Such a field can be created 



in the vicinity of a hard ferromagnetic nanostructure [19] , 
by applying a spatially uniform bias field that cancels the 
field of the nanostructure in a small region of space. In 
essence, the strong gradient field at one spin can substi- 
tute for strong anisotropic hyperfine couplings required 
for a purely molecular compass. This geometry provides 
a more significant anisotropy and thereby shows much 
larger directional sensitivity than does the conventional 
compass mechanism based only on anisotropic hyperfine 
couplings. Without requiring extra nuclear spins, the 
present model can work merely with two electron spins 
and thereby much simplifies quantum simulation of a 
chemical compass with, e.g. quantum dots and Nitrogen- 
vacancy centers in diamond. With more freedom to tune 
parameters in a better controllable environment, such 
kind of quantum simulations would be very helpful to 
understand the recombination process of radical pairs, in 
particular, whether and how quantum measurement and 
Zeno effect take place 0, [f| . 

Chemical compass mechanism. — Many chemical pro- 
cesses involve a radical pair intermediate, in which each 
radical has an unpaired electron coupled to an external 
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FIG. 1: (Color online) Left: A radical pair, coupled with the 
surrounding nuclear spins (black arrow), in a weak magnetic 
field B to be measured (yellow) and a strong magnetic gra- 
dient La (blue), due to e.g. a magnetic nanostructure. The 
outcome of a reaction depends on the direction of the weak 
field B. Right: The directions of B and the gradient field 
at the location of the acceptor La depicted in the molecular 
coordinate frame. 
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magnetic field and a few nuclei via the Hamiltonian [2C 



H 



k,j 



Sk ■ Xk 



(1) 



where 7 e = —g e fiB is the electron gyromagnetic ratio, 
Afej denote the hyperfine coupling tensors and Sk, Ik } 
are the electron and nuclear spin operators respectively. 
In our model, the magnetic field consists of two parts: 
Bk = B + Lfc, where the directional information about 
B is what one wants to infer from the radical pair reac- 
tion, and Lk is the local gradient field applied to each 
radical and is independent of B. The spin relaxation 
and decoherence times resulting from the factors other 
than hyperfine interactions are assumed to be consid- 
erably longer than the radical pair lifetime \3 L I 111 , to 
maximize sensitivity to weak magnetic fields [21] . In 
many photochemical processes, the radical pair is cre- 
ated in a spin-correlated electronic singlet state \S) = 
^7j(|t4-) — lit)) within the timescale of picoseconds. The 
nuclear spins start at thermal equilibrium, which under 
ambient conditions leads to an approximate density ma- 
trix as Pn(0) = &)j^j/dj, where dj is the dimension of 
the jth nuclear spin and Ij is the identity matrix. The 
Zeeman splitting from a magnetic field B as weak as the 
geomagnetic field is much smaller than the thermal en- 
ergy at ambient temperature. Nonetheless, the field can 
influence the non-equilibrium electron spin dynamics and 
thereby determine the ratio of the chemical product from 
the singlet or triplet recombination as long as the ther- 
malization time is longer than the reaction time. 

In experiments, one may measure different quantities 
that are dependent on the weak magnetic field B. Here 
we consider a simple first-order recombination reaction 
of the singlet radical pairs. We note that there is some 
controversy over how to describe the radical pair reac- 
tions (see e.g. 0, H, EH E3|)- Nevertheless, the con- 
ventional phenomenological density matrix approach [20| 
works well in most cases, in particular when the sin- 
glet and triplet recombination rates are the same (i.e. 
ks = kx = k) Q. We adopt this method and cal- 
culate the singlet yield as $s = J °° f(t)Ps(t)dt, where 
f(t) = ke~ kt is the radical reencounter probability distri- 
bution, and Psit) = (§|ps (t)|§) is the singlet fidelity for 
the electron spin state p s (t) at time t. The i ntegrat ion of 
$s was performed following the method in |25l |26| . 

Gradient enhancement of magnetic field sensitivity. — 
We starts from an optimally designed hyperfine compass 
model, one radical has strong and anisotropic hyperfine 
interactions, and the other radical has no hyperfine cou- 
plings |2]J . We arbitrarily choose to call the first radical 
the acceptor, A, and the second the donor, D, though 
nothing that follows depends on this designation. Ritz 
and coworkers proposed that the radical pair FADH -O^ - 
meets this criterion, and they further speculated that this 
radical pair may be responsible for the magnetoreception 



of European robins 2l| , see also 27| ■ Without loss of the 
essential physics, we take the hyperfine couplings (~ G) 



from FADH -Oj [28( for our calculations. 
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FIG. 2: (Color online) Magnetic field sensitivity of a chemical 
compass enhanced by a gradient field, (a) Singlet yield $s 
as a function of the angle 9 of the weak magnetic field B 
(B = 0.46 G) with different gradient field strengths on the 
acceptor, i.e. La — G (red, • • • o • • • ), 20 G (blue, ■■■<>■■■), 
40 G (green, • • • A • ■ • ), 80 G (purple, ■••*•••), while L D = 0. 
The recombination rate k = 0.5/is -1 . (b) Visibility V as a 
function of the radical pair lifetime r = 1/k. The direction of 
the gradient field La is set as 9a = 0. The same values of La 
are used as in (a). 

We define the molecular frame as the coordinate sys- 
tem, and the weak magnetic field B can be represented 
as B = B(sin 9 cos </>, sin 6 sin </>, cos 9) . The gradient field 
induces different local fields on two radicals. We as- 
sume that the gradient field on the acceptor radical is 
La = La{sui9a,Q,cos6a) while Ld = for the donor 
radical. The strength of the weak magnetic field to be de- 
tected is the same as the geomagnetic field, i.e. B — 0.46 
G. To demonstrate the basic idea, we first consider = 0, 
and then generalize to arbitrary (f>. 

In Fig. [5] (a) , we plot the singlet yield as a function of 
the angle 9 of the weak magnetic field B with different 
gradient field strengths L A =0G, 20 G, 40 G, 80 G on the 
acceptor. In the case of La = 0, the directionality comes 
only from hyperfine anisotropy. The gradient field clearly 
enhances the amplitude of the direction-dependent com- 
ponent of the magnetic field effect (MFE). To quantify 
the directional sensitivity, we use the magnetic visibility 
defined as Q 

V = (max<i>5 - min$s)/(max$s + min<i>s). (2) 

As the gradient field becomes larger, the sensitivity will 
increase and approach to a saturate best value. Fig. [5] 
(b) shows that for long radical pair lifetimes, the visibil- 
ity with the gradient field L A =40 G is almost twice the 
visibility without the gradient field. Usually, the radical 
pair lifetime should be very long (microseconds) to max- 
imize the effect of weak magnetic field [2l|, and hence 
performance, of the chemical compass [Fig. [2] (b)]. This 
requirement places a severe constraint on the chemistry; 
in typical radical pair reactions the lifetime is less than 
100 ns [20]. By increasing the overall magnitude of the 
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visibility, gradient-enhancement broadens the range of 
candidate reactions for a chemical compass. 

Liquid crystal experiment. — In a uniaxially oriented 
sample, the MFE is averaged over all values of the angle 
4>. Such a sample is prepared by, for instance, freezing 
the molecules in a nematic liquid crystal in the presence 
of a strong magnetic field |18j |. The ensemble-averaged 
MFE depends on only and is characterized by 



<M*)> 



1 

2^ 



271 



(3) 



It can be seen from Fig. [3] (a) that the enhancement of the 
sensitivity can still be observed with the average signal 
(<&s(0)) by choosing appropriate values of 9a- 

To induce the gradient field as above, one feasible way 
is to use magnetic nanostructures [l9j |. We model the 
nanocrystal as a uniformly magnetized sphere, in which 
case the external magnetic field is the same as that of a 
point dipole of magnetic moment m located at the center 
of the sphere [5^] . We denote the position relative to the 
center of the sphere by the vector r, and assume that 
both r and m lie along the z-axis. The magnetic field at 
r is 



B(r) 



2nr 3 



(4) 



where fio = in x 10~ 7 N-A~ 2 is the permeability of free 
space, the magnetic moment m — MgQ, with M the spe- 
cific magnetization, p is the material density, = ^R 3 
is the volume of the particle and R is its radius. The 
parameters for the typical magnetic material Fe,-;04 are 
M = 43 A-m 2 - kg" 1 , g = 5210 kg-m" 3 \V§\. For 
molecules with a separation rAD between two radicals 
a few nanometers 30], it is sufficient for a nanoparti- 
cle to induce a large local field imbalance (~ 10 G) on 
the donor and acceptor. For example, using a Fe3C>4 
nanoparticle with the radius R — 15 nm, it is possible to 
induce the local field difference as large as ~ 40 G be- 
tween the position = 35 nm and rrj = rA+rAD = 38.5 
nm (assuming rAD =3.5 nm). By generating an addi- 
tional homogenous field to compensate the field at the 
position re, we can effectively obtain the gradient field 
on the donor and acceptor molecule as La — 40 G and 
Ld = G respectively. 

To see whether the effect of the gradient field shown 
above can manifest with experimental imperfections, we 
take into account the fluctuations of La and Ld by mod- 
cling the fluctuation as the three-dimensional Gaussian 



distribution f(Ai 
with a a 



■ exp (- 



(27TCT 2 ) 3 / 2 

2 G and cr D = 0.1 G. Therefore, 
average of (<M#)) in Eq. © is 



) (i = A,D) 
the ensemble 



$s(«)= / ($s(0)> \A A ,A D f(A A )f(A D )dA A dA D (5) 

where ($5(0)) \a a ,a. d is the average singlet yield when 
the local fields on the acceptor and donor molecules are 
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FIG. 3: (Color online) Magnetic field sensitivity in a liquid 
crystal experiment, (a). Visibility of the average singlet yield 
(5>s(0)) as a function of the angle 6 a of the gradient field 
La = 80 G. The blue dashed curve represents the visibility 
without the gradient field, (b) Ensemble average of the singlet 
yield &s{8) in Eq. ([5| from a Monte Carlo simulation of 2x 10 4 
samples as a function of the angle 8 of the weak magnetic 
field B. The gradient field is La = 40 G with 6a =0 (purple, 
•••*••■), while Ld = 0. The fluctuations of the local fields 
Ld and La are characterized by the 3-dimensional Gaussian 
distributions with the variance a a = 2 G and <td = 0.1 G. For 
comparison, we plot the singlet yield with no gradient field 
(red, • ■ ■ o ■ ■ ■ ), and the one with the gradient field La = 40 G 
(6a =0) without fluctuations (blue, • ■ • A ■ ■ • ). In both panels 
the radical pair lifetime is 2(j,s and B — 0.46 G. 



La + Aa, Ld + Ad respectively. We have used Monte 
Carlo simulations to calculate the above ensemble aver- 
age in Eq.©. In Fig.[3](b), we see that the enhancement 
from the gradient field can still be observed. 

Probe spin correlations in a chemical compass. — Be- 
sides the significant enhancement of the directional sen- 
sitivity offered by gradient fields, we now examine how 
they can provide new insights into the quantum dynamics 
of radical pair reactions. For the present model chemical 
compass, if the gradient field on the acceptor La domi- 
nates over the hyperfine couplings and the weak magnetic 
field B, the singlet yield can be written as [2fJ 



(0) 



where the expectation value is calculated over the ini- 
tial state, and A = \uo)(uq\ — \ui){v,i\ (with {|uo)j l^i)} 
the eigen states of La ■ Sa), V = (U^AUd) with 
Ud = exp (i^j e tB ■ Sd) and the average taken over time 
weighted by f(t). By choosing La in the direction of x, 
y, and z, the corresponding operator A will be X, Y, Z 
respectively (which are the Pauli operators). Moreover, 
for each A, one can choose B also in the direction of x, 
y, and z such that the operators of V (as a linear com- 
bination of Pauli operators) are linear independent, see 
[26| . The singlet yields corresponding to these choices 
of La and B lead to nine independent equations, from 
which we can infer the spin correlations (M <g> N) for 
the radical pair state, where M,N = X, Y or Z. With 
these correlations, one may check whether the radical 
pair state violates Bell inequalities 3l|; or obtain lower 
entanglement bounds of the radical pair state, see Ref. 
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FIG. 4: (Color online) The gradient field as a tool to test the 
initial radical pair state. The singlet yield <I>s as a function of 
the angles (6, cj>) of the weak magnetic field B with B = 0.46 
G. The gradient field on the acceptor is La = 80 G, on the 
donor it is Lo — 0. The angle of the gradient field 6 a with 
respect to the z-axis is (upper), ? (middle) and ^ (lower). 
The patterns of the singlet yield over (6, <j>) from the initial 
singlet (left) and classically correlated state (right) are quite 
similar for 9a = 0, but are very different for 8a = j and ?. 
The radical pair lifetime is chosen as 2/is. 



32]. The above idea can be extended to monitor the 
dynamics of spin correlations suppose one can switch on 
gradient fields during the reaction. 

As an example, we show that gradient fields can dis- 
tinguish the singlet and the classically correlated initial 
state p c = (|tl) (til + lit) (ttl)/2- For systems where the 
radical pair lifetime is much longer than the decoherence 
time, the conventional hyperfine-mediated MFE does not 
strongly depend on the initial states and thus can not al- 
low one to achieve this goal, see e.g. 0). If the gradient 
field is along the z-axis, the singlet yields are quite simi- 
lar for the singlet and the classically correlated state [see 
Fig. 2] (upper)]. However, if we vary the direction of the 
gradient field La, then the visibility for the singlet state 
will be much larger than for the classically correlated ini- 
tial state [Fig. |4] (middle, lower)]. In particular, for the 
classically correlated state, the singlet yield is insensitive 
to the angles (9, </>) for 6 a = § while for the singlet state 
the angular sensitivity persists [Fig. [4] (lower)]. The dif- 
ference originates from the essential boundary between 
classical and quantum correlation (entanglement). The 
large gradient field can be viewed as a measurement of the 
acceptor spin along La- the singlet state demonstrates 



perfect anti-correlation of the spins for any direction of 
La, while for the classically correlated state this is true 
only in a certain direction of La (i.e. the z direction). 

Summary. — We have demonstrated that a gradient 
field can lead to a significant enhancement of the perfor- 
mance of a chemical compass. The gradient field also pro- 
vides us a powerful tool to investigate quantum dynamics 
of radical pair reactions in spin chemistry. In particular, 
it can distinguish whether the initial radical pair state is 
in the entangled singlet state or in the classically corre- 
lated state, even in the scenarios where such a goal could 
not be achieved before. These phenomena persist upon 
addition of partial orientational averaging and addition 
of realistic magnetic noise. The effects predicted here 
may be detectable in a hybrid system compass composed 
of magnetic nanoparticles and radical pairs in an oriented 
liquid crystalline host. Our work offers a simple method 
to design/simulate a biologically inspired weak magnetic 
field sensor based on the radical pair mechanism with a 
high sensitivity that may work at room temperature. 
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Appendix 

Calculation of singlet yield. — We adopt the method as in [l] to calculate the singlet yield. For the self- 
completeness, here we present a simple outline of this method. The Hamiltonian for the system (two electron spins, 
one of which is coupled with several surrounding nuclear spins ) is as follows 

h= J2 H k = -7e E ■ s k + E & • K ■ 4 ( 7 ) 

fc=l,2 fe k,j 

In our calculations of the main text, we have neglected the Zeeman interactions between the nuclear spins and the 
external magnetic field. We have included these interactions, and verified that the induced difference is very small 
(as the gyromagnetic ratio for a nucleus H and N is much smaller than 7 e ). The numbers of hyperfine couplings are 
take from Ref.[2|. We calculate the singlet yield as |3| 

/>oo 

$s = / f(t)P s (t)dt (8) 
Jo 

where f(t) — ke~ kt is the radical re-encounter probability distribution, and Ps(t) = (S\p(t)\S) is the singlet fidelity 
for the electron spin state p(t) at time t. Eq. ([5]) can be obtained from the conventional Haberkorn approach Q in the 
case that the singlet and triplet recombination rates are the same, i.e. k$ — kr — k. The singlet yield is calculated 
following the method in [![. We first write the singlet fidelity as 

P s (t) = Tr[e-^(p (g)^)e^(|5)(5|(g)I,)] (9) 

3 3 3 

= EEHO*)® ■ H(|5>(5|(g) 1,01m) . e -*"— )' (10) 

m n j > 3 j 

where po is the initial state of the radical pair, and the initial state of the nuclear spins at room temperature can be 
approximated as pb(0) = &)j^-j/dj, where dj is the dimension of the jth nuclear spin, and {|m)} and {|n)} denote 
the eigen states of the Hamiltonian H in Eq.(l). After some calculations, we have the singlet yield as 

$S = / ke- kt P S (t) = -EE PmnA nm —- (11) 

where d = ]Jdi , p m , n = (m\(p 0^- Ij)\n) and A n . m = (n\(\S){S\ ®j- lj)\m). 



Comparison between Haberkorn approach and quantum measurement master equation. — In the main 
text, we consider the radical pair reaction with the same singlet and triplet recombination rate, i.e. k$ = fcr = k. The 
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method we use to calculate the singlet yield is based on the Haberkorn approach 0] that describes the recombination 
of radical pairs 



^ = -i[H, p] - 'fiQsP + P Qs) - ^{QtP + P Qt) 



(12) 



where Qs and Qt are the projection operators for the singlet and triplet electronic states of the radical pair. There are 
alternative master equations based on quantum measurement that have been proposed to describe the recombination 
of radical pairs @, Under the condition ks — kx = k, these master equations [H, [|| can be written in the following 
form 



dp 



= -i[H, p] - (ks + k T )p + ksQrpQr + k T QspQs 



(13) 



From Eq. p^]) or Eq. (jT3")) . one can obtain the density matrix of the radical pair state pit) at time t and thus calculate 
the singlet yield as 



/>oo 

<S>s = k Tr [Q s p(t)} dt 
Jo 



(14) 



We compare the results of the singlet yield from these two approaches. For the simplicity of calculation, we take the 
three most significant hyperfine interactions in FAD IT, i.e. those for the nitrogens N5 and N10 and the proton H5 Q. 
It can be seem from Fig. [5] (for the long lifetime r = 1/k = 2ps) and Fig. [5] (for the short lifetime r = 1/k = 50ns) 
that the difference between the results from two approaches is very small (around 1 ~ 2%). 
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FIG. 5: (Color online) A comparison between the singlet yield from Haberkorn appraoch (HAB, red) and quantum measurement 
master equation (QM, blue). The radical pair life time is t = 1/k = 2ps. The other parameters are the same as Fig. 2 (a) in 
the main text. 



Probe spin correlations with gradient fields. — We assume that the gradient field on the acceptor La is much 
larger than the hyperfine couplings and the weak magnetic field B. To calculate the singlet yield, for simplicity, the 
Hamiltonian can be approximated as 0] 



H~- le (L A -S A + B-S D ) 



(15) 
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FIG. 6: (Color online) A comparison between the singlet yield from Haberkorn appraoch (HAB, red) and quantum measurement 
master equation (QM, blue). The radical pair life time is r = 1/k = 50ns. The other parameters are the same as Fig.2 (a) in 
the main text. 



We denote the density matrix of the initial radical pair state as p, and use the eigen states of La • Sa ({\uq}, \u\}}), 
namely La • Sa\uq) — La\uq) and La ■ Sa\ui) = — La\ui), as the spin basis of the acceptor. The initial state can then 
be written as 



P= \u m )A(u n \ ® Pd™ where p^' n = (u m \p\u n ) 

m.n— 0,1 



The singlet fidelity at time t is 



P s (t) = (S\ [\u )(u Q \ ® (U D f/$U* D j\ \S) + (S\ [\ Ul )( Ul \ <g> (UDp^ujj] \S) 

+ e a ^ tLA {S\ f|uo>(ui| ® (Uop'Bul)] \S) + e~ a ^ tLA (S\ \\ Ul ){u \ ® (U D p^uj,)\ \S) 



(16) 



(17) 



where Ud = exp (ij e tB ■ Sd)- If La is very large, the last two terms (second line) in the above equation oscillate very 
fast and make no effective contribution to the singlet yield due to time average, thus the singlet yield will be 



o 



f(t)-{(S\ \u )(u Q \®(U D p°°ul) \S) + (S\ |ui)(ui|®(tf D /#0£ \S)}dt 



2 7o 



f^-Trip^U^iu^UD+p^UlMiuolU^dt 



/(<)'%- 



I A® W D AU D 



)p]dt 



(18) 
(19) 

(20) 
(21) 



where A = |u )(moI — an d V = J °° dtf{t)U' D AUB- By choosing La in the direction of x, y, and z, the 

corresponding operator A will be X, Y, Z respectively (which are the Pauli operators). As an example, we assume 
that A = Z, it can be seen that if we choose B in the direction of x, y, and z respectively, the corresponding operators 
of V are cZ - sY, cZ + sX, Z, with c = J °° dtf(i) cos (j e tB) = k 2 /[k 2 + (-f e B) 2 ] and s = dtf(t) sin (^tB) = 
k~f e B/[k 2 + (-y e B) 2 }. These operators are linear independent. From the singlet yields corresponding to these choices 
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of La and B, we can have three independent equations, following which we can infer the spin correlations (Z ® X), 
{Z ®Y), (Z ® Z). In a similar way, we can choose A = X, Y and obtain the spin correlations (X <g> X), (X <g> Y), 
{X®Z), and {Y®X), {Y ®Y), {Y ®Z). 
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